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Abstract. The effects of the propagation of particles which have a finite life-time and an 
according broad distribution in their mass spectrum are discussed in the context of a transport 
descriptions. In the first part some example cases of mesonic modes in nuclear matter at finite 
densities and temperatures are presented. These equilibrium calculations illustrate the dynam¬ 
ical range of spectral distributions to be adequately covered by non-equilibrium description of 
the dynamics of two nuclei colliding at high energies. The second part addresses the problem 
of transport descriptions which properly account for the damping width of the particles. A 
systematic and general gradient approximation is presented in the form of diagrammatic rules 
which permit to derive a self-consistent transport scheme from the Kadanoff-Baym equation. 
The scheme is conserving and thermodynamically consistent provided the self-energies are ob¬ 
tained within the <I>-derivable two-particle irreducible (2PI) method of Baym. The merits, the 
limitations and partial cures of the limitations of this transport scheme are discussed in detail. 


1. Scope of high energy nuclear collisions 

Two atomic nuclei collide at high relative energies. During a time span of less than 10“^^ sec. 
a highly compressed and heated interaction zone is formed with matter densities p that exceed 
several times the nuclear saturation density po of 0.16 nucleon per cubic femto-meter (fm^) and 
temperatures T expressed in energy units beyond 100 MeV, cf. Fig. 1, below [1]. In daily 
life units this corresponds to p > 10^® g/cm^ and T > 10^^ K and thus to the highest matter 
densities and temperatures ever realized in laboratory experiments. Such matter conditions 
prevailed during the first few micro-seconds of the early universe, occur in super nova explosions 
and their remnants, the neutron stars at somewhat lower temperatures. This form of matter 
can only be investigated in the laboratory through nuclear collisions where one collects the 
reaction products emitted from the highly excited collision zone in sophisticated multi-particle 
detector arrangements. For a recent and comprehensive experimental overview see e.g. ref. [2]. 
In such experiments one distinguishes strongly interacting reaction products, such as hadrons, 
essentially nucleons and pions but also hadrons with other flavors like strangeness or charm 
from solely electro-magnetically interacting probes such as photons or lepton pairs (electrons- 
positrons or di-muons). While the former can only give a clear view back to the situation where 
the particles were essentially frozen out, i.e. the strong interactions cease, the latter are classified 
as penetrating probes. Once created these particles pass the collision zone essentially without 
further interaction. Such penetrating probes permit to look deeply into the interaction zone and 
are thus interesting as messengers from the compressed and heated matter. 
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Figure 1. Collision of two gold nuclei at a bombarding energy 
per nucleon of 160 GeV. Contour plots of the densities in units of po 
at two time steps, resulting from a three-fluid hydrodynamical model [1]. 


The observation of di-leptons has received much attention over the last decade due to the 
observed enhancement of di-lepton pairs in the invariant mass region below the vector-meson 
masses [3,4] and the observed suppression of J/T-mesons in central nuclear collisions compared 
to proton-proton collisions. With the recently completed HADES di-electron spectrometer GSI 
continues this research at lower bombarding energies. Considerable progress in mass resolution 
achieved by new pixel detectors for high energy muons has enabled new experimental results 
with unprecedented mass resolution [5,6]. 

The endeavor to investigate such nuclear collisions in experiments at GSI, CERN or at 
the RHIC-collider at the Brookhaven laboratory (BNL) strongly relies on adequate theoretical 
descriptions of the dynamics of such reactions. While a fully quantum mechanical or even field 
theoretical description is off reach, simplified approaches are considered that cover different 
dynamical ranges of applicability. With focus on the equation of state (EoS) of the interacting 
matter hydrodynamical models were investigated. Since a full instantaneous stopping to local 
equilibrium is not compatible with experimental observations, various kinds of multi-fiuid models 
were considered to relax these constraints and to allow for a partial transparency between the 
different fluid components [1], cf. the example given in Fig. 1. Even further generalizations 
treated the abundances of the different constituents by non-equilibrium rate equations, while 
only the kinetic degrees of freedom and the pressure were considered to be in local equilibrium [7]. 
These approaches permit to address the issues of phase transitions e.g. between the hadronic 
and the quark-gluon plasma (QGP) phase and the underlying EoS. 

At the microscopic level kinetic transport descriptions are in use. They describe the 
constituents of the matter as classical particles subjected to Newtonian forces and to stochastic 
collisions among the particles in terms of generalized Boltzmann equations which adopt mean 
fields and the statistics (Fermi-Dirac/ Bose-Einstein) of the particles. Thereby the collision 
term is essentially taken from binary-collision cross-sections. The relevant degrees of freedom 
may significantly change during the collision process. Due to the collisions among the incident 
nucleons the latter do not only change their momenta but rather also may get excited to baryon 
resonances, further more mesons and mesonic resonances can be created and at sufficiently high 
energy density the sub-hadronic degrees of freedom, the quarks and gluons, become liberated 
forming a QGP. 

In such collisions one is confronted with the fact that the dynamics has to include particles 
like the Delta isobar or e.g. rho-meson resonances which have life-times of less than 2 fm/c 
already in the vacuum. The equivalently damping widths thus exceed values of 100 MeV. Also 
the collisional damping rates deduced from presently used transport codes are of the same order, 
whereas typical mean temperatures range between 50 to 150 MeV depending on beam energy. 
Thus, the damping widths of most of the constituents in the system can by no means be treated 
as a perturbation. As a consequence the mass spectrum of the particles in the dense matter is 
no longer characterized by a sharp delta-function but rather acquires a width due to collisions 











and decays. One thus comes to a picture which unifies resonances that have already a width in 
vacuum due to decay modes with the ’’states” of particles in dense matter, that obtain a width 
due to collisions (collisional broadening). 

The theoretical concepts for a proper many body description in terms of a real time non¬ 
equilibrium field theory were already devised by Schwinger [8], Kadanoff and Baym [9], and 
Keldysh [10] in the early sixties, extensions to relativistic plasmas followed by Bezzerides and 
DuBois [11]. First investigations of the quantum effects on the Boltzmann collision term were 
given by Danielewicz [12], the principal conceptual problems on the level of quantum field 
theory were investigated by Landsmann [13], while applications which seriously include the 
finite width of the particles in transport descriptions were carried out only in recent times, 
e.g. refs. [12, 14-22]. For resonances, e.g. the Delta resonance, it was natural to consider 
broad mass distributions and ad hoc recipes have been invented to include this in transport 
simulation models. However, many of these recipes were not correct as they violated some basic 
principles like detailed balance (see discussion in ref. [14]), and the description of resonances 
in dense matter has to be improved [19-26]. A further conceptual problem in self-consistent 
resummation schemes concerns the proper renormalization of the (hidden) divergent loops in 
the scheme. In this question significant progress has recently been achieved [27-29]. 

In this contribution the consequences of the propagation of particles with short life-times 
are discussed. First we investigate at the example of nuclear matter at various temperatures 
a strongly interacting system in equilibrium. In particular the spectral properties of pionic 
modes in matter are explored. These pionic modes couple to vector mesons which through their 
electromagnetic decay into di-leptons (pairs of electron-positron or muons) are considered as 
special messengers of the dense matter phase. These example cases illustrate the properties of 
particles with broad damping width in dense matter environment. In the final part we discuss 
how particles with such broad mass-width can be described consistently within a transport 
theoretical picture. 

We are going to argue that the Kadanoff-Baym equations in the first gradient approximation 
together with the <I>-functional method of Baym [30] provide a proper frame for kinetic 
description of systems of particles with a broad mass-width. To this end, we discuss relevant 
problems concerning charge and energy-momentum conservation, thermodynamic consistency, 
memory effects in the collision term and the growth of entropy in specific cases [23]. 

2. Equilibrium properties of particles and resonances in dense nuclear matter 

In nuclear physics quite a broad variety of concepts and methods are in use depending on the 
particular question raised. This is due to the non-perturbative nature of the interaction. For 
nuclear structure commonly non-relativistic concepts are used, nowadays with so-called realistic 
two-body interactions, supplemented with a three-body correction. There the essential question 
lies in the proper treatment of the short range and tensor correlations [31, 32] with unifying 
prospects in the low-energy limit of a derived effective interaction [33]. In such approaches the 
exchanged bosons which furnish the forces only implicitly enter the dynamics through the static 
two-body potentials. 

For the here presented treatment with focus on high energy excitations the exchange bosons 
become dynamical. Therefore we employ an effective field theory picture with nucleons, and the 
Delta-isobar resonance as the main baryonic matter fields interacting via pion exchange for the 
long-range part of the interaction. This we supplement by a repulsive Migdal-type two-body 
interaction of zero range in order to provide proper saturation and other properties of nuclei. 
Interested in penetrating probes we extend the model to include the lightest vector mesons in 
the second part of this section. These vector mesons are discussed as important messengers 
about the dense nuclear environment produced in nuclear collisions, since they can be observed 
through their electromagnetic decay into electron-positron or muon pairs. 



2.1. Pions and haryonic resonances 

For the description of the nuclear matter we consider the nucleon and the A(1232)-resonance 
as the main baryonic degrees of freedom. The pion is the lightest strongly interacting boson. 
Thereby it is commonly considered as a Goldstone boson related to the approximate chiral 
symmetry which would exist for massless pions. In lowest approximation in the pion-energy the 
underlying symmetry leads to p-wave couplings to the baryonic pseudo-vector currents together 
with a minor s-wave term. The latter we ignore. Thus one has a standard 3-vertex which couples 
the baryonic current to the pion. For the energies involved we approximate the Dirac spinor 
structure of the baryons by the non-relativistic limit in no sea approximation (ignoring anti- 
baryons), while the kinematics is kept in the relativistic form. The explicit analytic expressions 
and further details can be found in ref. [34],see earlier work in this context [35,36]. To this 
end the pion self-energy is given by simple particle-hole and Delta-hole loops, while the baryon 
self-energies involve a single rainbow type loop which includes one pion line. 

The self-consistent coupled Dyson scheme includes all terms resulting from a RPA treatment 
of the one-pion exchange interaction between the baryons. The latter accounts for the long- 
range attraction among the nucleons. Besides this one needs a repulsive interaction of short 
range between the nucleons in order to stabilize the matter. This interaction is commonly 
approximated by a zero range form with a 4-baryon vertex, initially introduced by Migdal [37]. 
Iterating the corresponding nucleon-hole and Delta-hole loop insertions to all orders the retarded 
pion self-energy 
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can still be obtained in closed form [38,39]. Here and in the following the boson self-energies 
are denoted by H. The gik denote the Migdal parameters and F{q) is a pion-baryon form-factor 
that affects the high energy-momentum behavior of the pion. The particle-hole loops are given 
by 



Likewise the nucleon and A-resonance self-energies are given by 



The corresponding self-consistent results for the pion spectral function are shown in the Figs. 
2 and 3. In the spectral function of the pion we see the effect of the coupling to the different 
baryonic excitations. The on-shell pole broadens and is shifted towards lower energies due to 
the net-attraction of the interaction in this kinematical region. The pion strongly mixes with 
the AN~^ component. The latter leads to the shoulders around the peak. This component 
appears in the time-like region for low momenta, while it traverses to the space-like region for 
higher momenta. At low energies an entirely space-like pion component appears due to NN~^ 
excitations. From the kinematics it is clear that this component stays in the space-like region 
for all momenta. Such space-like pion modes have to be interpreted as a scattering process 
between the baryons in the matter mediated by the exchange of pions. Due to this component, 
the pion spectral function has non-vanishing strength at all energies, such that in all processes 
where pions are involved all thresholds disappear. The short-range Migdal repulsion prevents 
pion condensation to occur at nuclear saturation density. The form of the spectral function 



turns out to be sensitive to values of the Migdal parameters and the pion baryon form-factor 
F{q) employed in the model. Due to this form-factor, which is fitted to the vriV-scattering data 
the A-hole component appears only as a slight shoulder above the peak. For the resolution used 
in the calculations the well known zero sound component at zero T is invisible. 

For higher temperatures all structures are broadened. It is clear that such struc¬ 
tures can no longer be described by a quasi-particle picture. Rather one has to take 
due account of the fully dressed spectral distributions of the particles in the medium. 



Figure 2. Contour plot of the loga¬ 
rithm of the pion spectral function A.^ at 
T=0 and 120 MeV (upper/lower panel) 
and p = po- The line spacing covers half a 
decade. 
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Figure 3. Pion spectral function 
at T=0 and 120 MeV (upper/lower panel) 
and p = pq for two different momenta. 


2.2. Vector-mesons coupled to the in-medium pion cloud 

In view of the strong modifications of the pion spectral function through the dense nuclear 
medium it is interesting to study its influence on other modes that strongly couple to the pion, 
like the vector mesons, i.e. the p- and w-mesons. The vacuum widths of these resonances are 
generated through the decay into two pions in the case of the p-meson, respectively three pions 
for the cu-meson. For the w-meson we choose the indirect decay via the p-meson (the so called 
Gell-Mann, Sharp and Wagner (GSW) process [40]) which is the dominant channel in vacuum. 
The coupling Lagrangian chosen in accordance with the low energy limit of QGD [41-43] gives 
rise to the following vector-meson polarization tensors (self-energies) 

nr (9) = (2) 
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with four-momentum q = (q'Ojq)- The two vector-meson coupling vertices are adjusted to the 
vrTT-scattering phase shifts together with the electromagnetic form-factor of the pion and to the 
vacuum decay width of the a;-meson, respectively. For the numerical results we neglected the 
coupling of the p-w-loop back to the pion as the vector mesons cause only a minor perturbation 









to the TrA^A-system. The correlation between p- and w-meson modes results from the cj-meson 
self-energy © and the reverse process encoded in the coupling of the p-meson to the w-vr- 
loop which both are included self-consistently. Since the model omits higher lying degrees 
of freedom the in-medium changes of the real parts of the vector-meson self-energies are less 
precisely determined within the model space. Therefore we drop them during the self-consistent 
iterations, while restoring the normalization of the spectral functions at each step. A posteriori, 
with the parameters used we checked for the thus ignored real parts of the loops using a counter¬ 
term scheme which through vector dominance fixes the photon mass (i.e. at = 0) to be zero 
with residue equal one in vacuum. It gave values of Ren/(2go) which are by far less than the 
corresponding widths of the vector mesons and therefore negligible. Keeping this subtraction 
scheme also at finite T, i.e. ignoring contributions from hidden divergences [27] and tadpoles, 
which solely lead to an additional T-dependent mass shift, we still obtained insignificant changes 
in this quantity. 

The resulting polarization tensors 0 have to be four-transversal because of current 
conservation. This is provided by a strategic projection method which uses two moments of 
the spacial part of the respective polarization tensors in the self-consistent loop approximation 
and constructs the respective four-transverse tensors. For details we refer to refs. [34,44]. 

In the following we explicitly analyze the influence of the various components of the pion 
spectral function on the vector-mesons and thus on the di-lepton spectra. Formally we do this 
by splitting the spectral function of the vector-meson into the various components related to 
the different processes feeding into this vector-meson channel. Thus, decomposing the total 
damping width into partial widths Fv,tot(Q) = the di-lepton yield can be brought 

into a Breit-Wigner like form with partial in- and out-widths 
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Figure 4. e'’'e“-rate from the decay of the 
cu-meson at T=120 MeV, p = pQ divided 
into different contributions: 
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Figure 5. e"’“e“-rate from the decay of the 
p-meson at T=120 MeV, p = po divided into 
different contributions: 
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(suppressing the tensor structure of spectral function and polarization tensor which leads to a 
degeneracy factor 3 for vector particles). Here denotes the vector meson mass, /t(qo) is 


















the thermal Bose-Einstein distribution and r^g+e- = “Iiiinve+e-/(2(i'o) is the di-lepton decay 
width of vector-meson v. In the vector dominance picture the vector mesons directly couple to 
a virtual photon which then decays into the respective lepton pairs, where for massless leptons 
Ev e+g- oc a‘^m^j{2qoq^) with fine-structure constant a. 

We start the discussion of the effects on the cu-meson. In the medium there are three major 
processes contributing to its damping width, illustrated by perturbative “time-flow” diagrams 
where the time is running from left to right and vertical lines denote virtual space-like propaga¬ 
tors, which mediate two-body interactions: 


P UJ P LO 



TT N N 


The subsequent decay of the tu-meson into the virtual time-like photon and its final decay into 
the lepton pair is not illustrated. In the self-consistent calculation with finite width spectral 
functions all these processes are included automatically by one self-energy diagram. For the 
hrst fusion-type process (F), p7r —> cj —>■ e"'“e“, the cj-meson is formed by the fusion of a p- 
meson with a time-like in-medium pion. Its inverse exists already in vacuum and determines the 
vacuum decay width of the cj-meson. The second process (D), p —ttcj ^ 7re'''e“, corresponds 
to a p-Dalitz-decay via an intermediate tu-meson. In the self-consistent calculations, both above 
mentioned processes just differ in the sign of the pion energy in the vrp-loop of the cu-self-energy 
©• The process (S) in (jH) corresponds to the scattering pN—> wN mediated by a virtual, 
i.e. space-like pion exchange. In view of the pion modes at zero temperature (cf. Fig. 2) we 
isolate this space-like component by a cut on the far space-like region with pion loop momenta 
with IpI > 2|po| in ©. At T = 120 MeV this separation is less evident in view of the broad 
structure of the pion spectral function (Fig. 3). Thus the different components of the processes 
displayed in Fig. 4 somewhat depend on this cut. In addition to the fusion width (F), which 
constitutes just a temperature dependent modification of the vacuum width (e.g. accounted for 
by Schneider and Weise [45]) a genuine in-medium process, namely the scattering process (S), 
contributes with comparable strength at the nominal resonance position. The “Daliz”-decay of 
the p-meson (D) dominates the low mass region. 

In summary of this analysis: a major portion to the u) spectrum results from processes (S) 
which are not accounted for in the simple on-shell treatment. These contributions, however, 
sensitively depend on the in-medium properties of the virtual pion cloud, which certainly needs 
further clarifying investigations before quantitative conclusions can be drawn. 

The same type of processes as in (@]) also occur for the p-meson just interchanging p with 
UJ, listed in the p-meson decomposition given in Fig. 5 as (F), (D) and (S^^). However, these 
w-induced components are less important compared to the coupling to the two-pion channels 
given by 

vr P TT P 
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At invariant masses above 300 MeV, the 7r'’'7r“-annihilation process (A) is clearly dominant. 
Processes (B) and (S"^) are not present in any on-shell treatment of the pion as they solely 
arise from genuine off-shell components of the pion spectral function. In this respect (B) 
can be interpreted as a bremsstrahlung process radiated off a pion scattered in the medium. 


while process (S) corresponds to inelastic ttN —> pN scatterings mediated via virtual pion 
exchange. The latter two components, which emerge completely consistently within the model, 
only contribute to the very low-mass part of the invariant mass spectrum. 

3. Conserving non-equilibrium dynamics of particles and resonances 

The equilibrium results discussed above can immediately be implemented into a hydrodynamical 
scheme which relies on the assumption of local equilibrium. This concerns the EoS derived from 
the corresponding thermodynamic potential, as well as the spectral properties of the particles 
and resonances as the electromagnetic decay rates. As penetrating probes the latter can simply 
be integrated over the hydrodynamic space-time evolution, while strongly interacting probes 
require an appropriate freeze-out concept. 

On the microscopic level the Kadanoff-Baym (KB) equations [9] provide the appropriate 
framework for a self-consistent dynamical treatment. However, only very recently first attempts 
to solve these equations for very simple scalar models were successfully under taken [46, 47] 
in admittedly very simplified cases of simple scalar fields in spatially homogeneous systems. 
For the complex situations in nuclear dynamics substantial progress in computing power is 
required before meaningful calculation with the KB-equations can be performed. Therefore most 
application are dealing with transport approaches where a Boltzmann-like transport equation is 
solved for the involved microscopic degrees of freedom. The numerical solution is mostly done 
by Monte Carlo simulation methods. 

In this context one of the key issues was to generalize the concepts such that they can also 
properly deal with resonances or more generally with particles which have non-trivial i.e. non 
quasi-particle like spectral functions. The corresponding concepts were already laid out in the 
book of Kadanoff and Baym [9]. It involves a systematic gradient approximation of the Wigner 
transformed KB-equations, avoiding the step towards the quasi-particle approximation. Thereby 
the question of thermodynamic consistence and conservation laws is of vital importance. Is such 
a transport scheme thermodynamically consistent and are the conservation laws associated with 
the symmetries of the underlying theory fulfilled? 


3.1. ^-derivable scheme and exact conservation laws 

For self-consistent Dyson-resummation schemes and thus for the KB-equations the question of 
thermodynamic consistence and conservation laws was settled by Baym [30]. He showed that 
such schemes are conserving, if and only if all self-energies result from the functional variation 
of a functional, called which is a functional of the self-consistent classical fields (one-point 
functions) and two-point propagators, see also [23]. The functional method initially introduced 
by Luttinger and Ward [48] and later reformulated in terms of path integrals became known as 
the CJT-formalism [49] in field theory according to the names of the authors. The 4>-functional 
is given by two-particle irreducible (2PI) closed diagrams in accordance with the Lagrangian of 
the system expressed in terms of the self-consistent propagators and classical fields, while the 
vertices remain the bare ones. 

With focus on the Kadanof-Baym equations we discard the dynamics of the classical field for 
the further considerations and concentrate on the self-consistent dynamics of the propagators 
G{x,y). In a 4>-derivable scheme the self-energies are generated from the functional <I>{G, A} 
through the following functional variation^, cf. [23,30] 


-iF(x,y) = T 


(5i^{G,A} 

SiG{y,x) 


for real fields 
for complex fields 


where 
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(7) 



(upper/lower signs refer to fermions/bosons).Here and in the following labels for the different 
species in the system and internal quantum numbers are suppressed. Taking all 2PI diagrams for 
provides the exact self-energy and two-point propagator of the theory. The space dependent 
coupling factor A(x), whose physical value is unity, is introduced to permit further variations, 
cf. Eq. (jHl) below. 

The virtue of the <l>-derivable scheme is, that the diagrammatic series 0 of can be 
truncated at any level, without spoiling the conserving properties of the self-consistent scheme 
for the Dyson-Schwinger or KB equations, see Eq. m below. The so constructed self¬ 
energies lead to a coupling between the different species in the system, which obey detailed 
balance generalizing the special recipes for broad resonances given in ref. [50]. It has been 
shown [30], see also [23], that such a self-consistent scheme is thermodynamically consistent and 
exactly conserving at the expectation value level of the conserved currents. 

For local couplings the interaction energy density is given by a functional variation with 
respect to the interaction strength X{x) 
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The single-particle potential energy density is defined as 
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fP°t(x) = ^ f d'‘7/[S(x,y)(=Fi)G(7/,x)-h (Ti)G(x,y)S(y,x)] (9) 

J C- 

here written for complex fields. Both, and £’p°*, enter the energy-momentum tensor (I12j) . 

The <I>-derivable scheme ® implies exact conservation laws for the Noether currents and the 
energy-momentum tensor [23] given by 


d^.P{X) = 0, = 0 with (10) 

^ I ^eu^(^i)F(X,p), (11) 

= / ^^V(Ti)i"(^,p) + (f“*(^) - £^^\X)) . (12) 

Here 

F{X,p) = ^iG-^{X,p) (alias TiG<(X,p)) (13) 

denotes the distribution function in four-phase-space given by the four-Wigner transformation 
of the two-point propagators G ^{xi,X 2 ), while e stands for any of the conserved charges and 
the kinematical factor unrelated to the four-velocity is define in m below. 


3.2. Kadanoff-Baym equations and complete gradient approximation 

We assume the reader to be familiar with the real-time formulation of non-equilibrium field 
theory on the so called closed time contour, ubiquitously used in this book. Since we will deal 
with general multi-point functions we use the more convenient contour-vertex notation of 

refs. [23,51] rather than the more clumsy “< >” notation. In the following we follow the lines 
given in ref. [52]. The set of coupled KB equations on the time contour in “—h” notation^ reads 

(G-i(-i9i) - G^\-id2)) G'-+(l, 2) = ^ d3 (S(l-, 3)G(3,2+) - ^(1“,3)S(3,2+)) 

^ The numbers 1, 2 and 3 provide short hand notation for space-time coordinates xi,X 2 and xs, respectively, 
inclnding internal quantum numbers. With superscript, like 1“ and 2"*", assigned to them, they denote 
contour coordinates with — and + specifying the placement on the time or anti-time ordered branch. For 
an arbitrary two-point function / its decompositions into the two branches of the contour are denoted as 
/*^*(1,2) = /(1*,2*) with k,l £ The match to the notation used, e.g., in refs. [9,12] is given by 

/■+ = /<; /+■ = />; /■■ = r; /++ = /“• 
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and likewise for G~^ . Here 


Go^ip) = 


2 2 
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for relativistic bosons 

for non-rel. fermions or bosons. 
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is the inverse free Green function in momentum representation. The upper/lower sign factors 
refer to fermions or bosons, respectively, Go and G correspondingly denote the free and full Green 
functions. The driving term on the r.h.s. of Eq. m, summarized by G, is a functional of the 
Green functions through the self-energies S contour folded with G. The real-time integration 
contour is denoted by C. 

As usual the step towards transport equations is provided by introducing the four-dimensional 
Wigner transforms for any two-point function / through 


fix,y) = 


d^P 


(2vr) 


or its inverse relation. The KB Eq. m then transforms to 
vf^d^F{X,p) = G-+{X,p-,{G}) 


d 


with F(X,p) = ^iG ~^(X,p) and = t/^Gq ^(p). 


(16) 

(17) 

(18) 


Here the r.h.s. of the KB equation is also expressed in terms of the Wigner transforms of all 
Green functions through |T6|). The final step is to expand the complicated r.h.s. of Eq. (dzi) 
to the first-order in space-time gradients, cf. [Appendix A| or [52] for more details. Then the 
local part of this r.h.s., G^^^, defines the collision term. It entirely consists of non-gradient 
terms, where all the different mean positions (x* + Xj)/2 occurring in the various Wigner-Green 
functions are replaced by the externally given mean position X of the l.h.s., i.e. X = (xi-\-X2)/2 
and the diagrams are to be evaluated as in momentum representation. The corrections for the 
displacement to the true coordinates of each Green function are then accounted for to first order 
in space-time gradients. Here we simply abbreviate the gradient terms by a <0 operator defined 
Appendix Aj acting on the local diagram expression 


m 


where 


vf^d,iTi)G-+{X,p) = (1 + iO) {G(7 +)(X,p)} 


G(7+)(X,p) = G-+(X,p;{G(ioc)}), 


(19) 


( 20 ) 


is a functional of the local Green functions G(io(,) = G{X,p). Here and for all further 
considerations below, both. Green functions G{X,p) and self-energies S(A, p), whenever quoted 
in their Wigner function form, are taken in local approximation, i.e. with X given by the 
external coordinate and T,{X,p) void of any gradient correction terms. The gradient terms (seem 
[Appendix A ) have the form that in the diagrams defining Gqoc) pairs of Green functions are 
replaced by their space-time and momentum derivatives, respectively, just leading to equations 
linear in space-time gradients. Naturally the result of the 0-operation depends on the explicit 
form, i.e. diagrammatic structure, of the functional on which it operates. For the non-gradient 
term in Eq. m, which defines the local collision term 


G(7+) = ^J:-\X,p)akiG^+{X,p) - {T)G-^{X,p)akiJ^^+{X,p) 

diagram 

= TiE-+(X,p)iG+-(X,p) - (^)iG-+(X,p)iS+-(X,p), 

value \ 


( 21 ) 

( 22 ) 


gam 


loss 















we give both, the diagram expression (EU) and the normally quoted value expression ((22l). The 
latter simplifies due to a cancellation of terms which however survive for the order sensitive 
gradient operation. Here aik = = diag(l,—1) defines the “contour metric”, which accounts 

for the integration sense, and summation over the contour labels k, I € , +} is implied. 

The above quantum kinetic equation (fT9|) has to be supplemented by a local Dyson equation 
for the retarded Green function [9] 

{G^{X,p))-^ = (Go«(p))“' - S«(X,p), (23) 

which together with Eq. m provides the simultaneous solution to G~^ . The full retarded 
Green function G^ depends on the retarded self-energy = S — S ^ again 

in local approximation. Gq is the free retarded Green function. Please note that equation (P|) is 
just algebraic although it is obtained within the frame of the first-order gradient approximation. 

In most presentations of the gradient approximation to the KB equations, Eq. (IM is 
rewritten such that the gradient terms are subdivided into different physical parts: on the 
one hand Poisson bracket terms describing drag- and back-flow effects which are pulled to the 
left side of the equation and in cases where the self-energy explicitly contains internal vertices 
an additional memory collision term 

»'‘flpF(X,j.) + {ReE«,±iO-+} + {±iE-+.ReO«} =C-+,(X,p) + C-„t„,(X.p). (24) 

"th = S-„G,(X,p)^HG'+(X,p)-G-‘(.Y,p)a„E'+.„yY,p) 

0-+(X.p)Ej;.„,(X.p),(25) 

where is only non-zero if internal vertices exist in S. 

Depending on the questions raised either the here presented separation of the gradient terms 
or the compact formulation (HU) is more advantageous. For the derivation of the conservation 
laws only a unified treatment of both, the Poisson brackets and memory collision terms, reveals 
the symmetry among these terms, which then displays the necessary cancellation of certain 
contributions such that the conservation laws emerge. Starting form Eq. (fT9|) and using the 
diagrammatic gradient rules given in [Appendix A| one could indeed prove that a) the Noether 
currents m and b) the local approximation of the energy-momentum tensor (II 2|l are exactly 
conserved. In particular for the latter derivation the diagram rules (1X2)) to (IXtI) where an 
essential help to isolate the proper counting of diagrams with terms proportional to the number 
of vertices contributing to £’int in m and those proportional to the number of Green functions 
contributing to Spot- 

3.3. Physical interpretation of terms 

For the physical interpretation it is advantageous to introduce the spectral function A{X,p) = 
—21m.G^{X,p) determined by the retarded equation (OHl) . as well as the four-phase-space 
occupation functions f{X,p), by means of F{X,p) = f{X,p)A{X,p). In thermal equilibrium 
f{X,p) becomes a Fermi-Dirac or Bose-Einstein distribution in the particles energy pQ. The 
evolution of F is governed by transport Eq. (OH) or equivalently by (fTU)). Together with the 
retarded equation (EHI) this defines a generalized quantum transport scheme which is void of 
the usual quasi-particle assumption. The space-time evolution is completely determined by the 
initial values of the Green functions at time zero for each space point. Thus the evolution is 
“'MarkoviaF'', since the memory part of the collision term is kept up to first-order gradient 
terms only. Within its validity range this transport scheme is capable to describe slow space- 
time evolutions of particles with broad damping width, such as resonances, within a transport 
dynamics, now necessarily formulated in the four-dimensional phase-space. 




Coming from the usual on-shell Boltzmann or Boltzmann-Uehling-Uhlenbeck collision term, 
each occurring three-momentum distribution function together with its momentum integration 
is simply to be replaced by its four-momentum analog, thus 


f{X,p) 


d^p 


F{X,p) 


d^p 


f{X,p)A{X,p) 


d'^p 


(26) 


Alongside the normally occurring two-body cross-sections have to be replaced by the 
corresponding T-matrix expressions providing the proper off-shell extensions. For genuine 
momentum dependent T-matrices the collision term has a hnite virial due to the interactions 
at finite distances and therefore the collision term contributes to the conservation laws in a 
non-trivial fashion. Within field theory model applications with local couplings one simply has 
to evaluate the corresponding self-energies for the collision term. In this case the local collision 
term drops out of the conservation laws. With <I>-derivable self-energies the collision term assures 
detailed balance. 

More subtle are the first order gradient terms given by the two Poisson brackets and the 
memory collision term in Eq. (El. All three terms contribute to the conservation laws. Thereby 
the first Poisson bracket furnishes the so-called drag-flow. In the quasi-particle limit it accounts 
for the dressing of the particles by the dragged matter cloud as to form a quasi-particle with 
a non-trivial dispersion relation with a corresponding in-medium group velocity that can be 
expressed by an effective mass. This change in flow is just compensated by the second Poisson 
bracket through the polarization of the medium. The latter thus forms a back-flow component. 
Only the coherent play of both Poisson brackets restores the conserved Noether currents and 
thus recovers e.g. Galilei (Lorentz) invariance. 

Since the first Poisson bracket involves space-time and momentum gradients directly acting 
on the distribution function F this term has an easy classical interpretation where the motion 
of the corresponding particle is subjected to a force which generally is momentum dependent. 
A generalization of this concept to the four-momentum picture is straight forward, since it 
just amounts to establish the corresponding characteristic curves of the homogeneous hrst-order 
differential equation. For the second Poisson bracket term on the other hand the derivatives 
of the distribution function appear only implicitly through the self-energy with the result that 
they affect momenta other than the momentum externally entering the transport equation. This 
has to be such, since the discussed term describes the reaction of the surrounding matter on 
the particle moving through the matter. However this term escapes an immediate description 
in terms of test particles, such that a simulation algorithm could not yet be established for the 
exact quantum kinetic equation (EH). 

Guided by equilibrium relations Botermanns and Malfliet [15] suggested a simplihcation of 
this second Poisson term, cf. [21,22,24] 

{-iS-+,ReG^} ^ {/(A:,p)r(A,p),ReG^} with r(A:,p) =-2ImS^, (27) 

BM 


formally valid up to second order gradient terms. Here the distribution function f{X,p) directly 
appears, while F is the damping width. The advantage of this substitution is that now the 
Poisson-bracket derivatives directly act on the distribution function / and the term amends a 
test-particle simulation [21,22]. The price to be payed is that then the conservation laws are 
slightly modified, since instead of the spectral function A rather the entropy-spectral function 
A, = iAP^, as introduced in ref. [53], enters the conserved current expression 

4j^(X) = I ■0^ev>^f{X,p)As{X,p){X,p). 


( 28 ) 






The BM-substitution accounts for part of the back-flow. In the quasi-particle limit both spectral 
functions converge to the same d-function at the quasi-particle energy [24,53]. A further merit of 
the BM-substitution is that for certain collision terms an entropy current can be derived which 
fulfills an exact H-theorem, for details see ref. [24]. 

4. Summary and prospects 

The initially presented nuclear physics examples show that in the dense nuclear environment 
the particles acquire strongly broadened spectral functions. They reach damping widths which 
are of the same order as the systems temperature or even larger. Such spectral functions escape 
the conventional treatment within the quasi-particle limit. 

For the application to nuclear-collisions dynamics problems the Kadanoff-Baym equations 
within a certain truncation level would certainly be the ultimate description level. However 
the required numerical afford presently limits such application to simple space-homogeneous 
systems and prevents realization in full three space dimension, time and four momentum space! 
Therefore we review the prospects of generalized transport equations derived in first-order 
gradient approximation from the KB-equations. Avoiding the quasi-particle limit the resulting 
set of dynamical equations consists of a generalized transport equation and a retarded equation 
for each type of particle in the system. 

We conclude that the resulting equations posses remarkable properties if derived in a so called 
<I>-derivable scheme. There the driving terms, the self-energies, are generated from a functional of 
closed diagram, the <I>-functional, of the non-equilibrium Green functions that can be truncated 
at any level. Thereby the transport equation consists of a local collision term which obeys 
detailed balance and first order gradient terms. In their exact form the first order gradient terms 
guarantee exactly conserved Noether currents and a conserved energy-momentum tensor (the 
latter provided that local couplings are used). In thermal equilibrium they furnish a non-trivial 
equation of state in accordance with the truncation level of <h. The equations of motion further 
preserve the retarded properties of the real-time components of the contour Green functions. 
Thus, the dynamical equations provide a generic scheme for a proper treatment of the motion 
of particles and resonances with non-trivial, i.e. finite mass-width spectral functions. Gertain 
gradient terms escape an immediate numerical simulation in terms of classical test particles. 
However first numerical applications [46] with the Botermanns-Malfliet substitution [15,22,24] 
look promising. 

Initially we aimed at a proper transport description for broad resonances. However the 
presented schemes with dynamical spectral functions offers even more. Already at the lowest 
order self-energy level (genuine two-point function with no further internal vertices), which leads 
to a collision term and therefore to a finite damping width, the scheme provides the following 
advantages: 

• the scheme includes higher order processes such as Bremsstrahlung which would be absent 
at the very same diagram level in the quasi-particle approximation; 

• it leads to completely regular expressions at any order of the self-energy; in order to avoid 
double counting and the inclusion of redundant processes which are generated from lower 
order processes through the iteration of the equations of motion, higher order self-energy 
terms have just simply to be derived from a 2PI generating functional; 

• therefore it avoids serious conceptual problems that arise in the use of the quasi-particle 
approximation for higher order processes (multi-particle scattering). The latter lead to 
mathematically ill defined expressions (squares of delta-functions or principle values) for 
the collision term which precisely arise whenever an intermediate propagator kinematically 
acquires on-shell conditions. The latter processes are precisely the reducible ones, as they 
can be generated through an iteration of lower order processes. Within the quasi-particle 



picture there is no clear cut way to avoid such complications and mathematical pathologies. 
Only partial resummation schemes such as the Dyson resummation with finite damping 
width spectral functions can cure this problem. 

The partial resummations implied by the self-consistent Dyson or KB scheme lead to violations 
of Ward-Takahashi identities for higher order correlation functions. In the context of vector 
mesons this affects the current conservation of the polarization tensor. As we briefly discussed, 
special projection strategies [34, 44] can be used to recover the four transversality. Also in 
cases of spontaneously broken symmetries the Nambu-Goldstone theorem is violated as the self- 
consistent Goldstone-boson masses are non zero. Recently the authors [54] suggested a correction 
term to the <h-functional which recovers massless Goldstone bosons in the broken phase. 

In conclusion, the generalized transport equations together with the <h-derivable concept 
represent a generic dynamical scheme to deal with particles which attain broad spectral 
functions. In addition to the merits of this approach there are prospects that some of the 
deficiencies can be cured by adapted corrections without leaving this conceptual frame. E.g. 
with applications to Quantum Chromodynamics (QCD) in mind this gives hope that at a proper 
time also non-trivial self-energies can be treated for the case of gauge fields. 
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Appendix A. Diagram rnles for the gradient approximation 

Let M(l, 2) be any two-point function with complicated internal structure. We are looking for its 
Wigner function M{X,p) with X = \{xi +X 2 ) to first-order gradient approximation. The zero- 
order term is just given by evaluating M(l, 2) with the Wigner functions of all Green functions 
taken at the same space-time point X = (xi +X2)/2 and the momentum integrations being done 
as in the momentum representation of a homogeneous system. To access the gradient terms 
related to any Green function G{i,j) involved in M(l,2), its Wigner function G{^{xi + Xj),p) 
is to be Taylor expanded with respect to the space-time reference point X = {xi + X2)l2, i.e. 


G(^^,p)^G{X,p) + \ 




d 


\ dXt^ 


G{X,p). 


(A.l) 


Both, dxG{X,p) and the accompanying factors {xi — xi) and {xj — X 2 ), can be taken as special 
two-point functions, to which special diagrams can be assigned 


= \{di + dj) G{i,j) dxG{X,p), (A.2) 

i j 

-i{xi-Xj) —>-(27r)^^(5(p) (A.3) 

with the corresponding Wigner functions on the right hand side. Then the gradient terms of 
any complicated two-point function (given in different notations) can graphically be represented 
by the following two diagrams on the r.h.s. 


0{M(1,2)} = 0 

1 


2 




1 2 


3^4 3^4 


i 

\ 

M' 

+ 

M' 

1 


2 1 



(A.4) 




Here the diamond operator 0 formally defines the spatial gradient approximation of the two- 
point function M to its right with respect to the two external points (1, 2) displayed by full dots. 
The diagrammatic rules are then the following. For any G{3, 4) in M, take the spatial derivative 
dxG{X,p) (double line) and construct the two diagrams, where external point 1 is linked to 3 
by an oriented dashed line, and where point 2 is linked to 4, respectively. Interchange of these 
links provides the same result. Here M' is a four-point function generated by opening M(l,2) 
with respect to any propagator G(3,4), i.e. 


M'(l,2; 


3, 


4)=T 


(5M(1,2) 

(5iG'(4,3)' 


(A.5) 


The diagrams in Eq. (IA.4I1 are to be evaluated in the local approximation, i.e. with all Wigner 
Green functions taken at the externally given space-time point X. The dashed line (jQ adds 
a new loop to the diagram, which, if integrated, leads to momentum derivatives of the Green 
functions involved in that loop. Both, double and dashed lines have four-vector properties, and 
the rule implies four-scalar products between them. From (ESI) follows that 


1 3 





(A.6) 


In momentum representation these rules correspond to the partial integration. They also imply 
0{Af(l,2)} = 0, if M contains no internal vertices. The convolution of two two-point functions 



= {A{X,p), B{X,p)} + A{X,p)0{B{X,p)} + 0{A{X,p)}B{X, P), (A.8) 

where C(l,2) = J d3A(l, 3)H(3, 2). 


Besides the Poisson bracket expression {A,B} it leads to further gradients within each of the 
two functions. Applied to the r.h.s. of Eq. cni), this rule indeed provides the decomposition 
into Poisson bracket and memory terms of the quantum kinetic equation 1211 ). 
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